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HIGHLIGHTS 


►  Advanced  SOC  and  SOH  estimation  of  a  lithium-ion  battery  by  a  Dual  Kalman  filter  in  combination  with  support  vector  machine. 

►  Feasibility  proof  and  validation  by  measured  drive  cycle  data. 

►  Estimation  accuracy  far  better  than  the  reported  ones  (SOC  error  <  0.5%). 

►  Real-time  diagnosis  by  management  system  conceivable. 
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Two  novel  methods  to  estimate  the  state  of  charge  (SOC)  and  state  of  health  (SOH)  of  a  lithium-ion 
battery  are  presented.  Based  on  a  detailed  deduction,  a  dual  filter  consisting  of  an  interaction  of 
a  standard  Kalman  filter  and  an  Unscented  Kalman  filter  is  proposed  in  order  to  predict  internal  battery 
states.  In  addition,  a  support  vector  machine  (SVM)  algorithm  is  implemented  and  coupled  with  the  dual 
filter.  Both  methods  are  verified  and  validated  by  cell  measurements  in  form  of  cycle  profiles  as  well  as 
storage  and  cycle  ageing  tests.  A  SOC  estimation  error  below  \%  and  accurate  resistance  determination 
are  presented. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  automotive  applications,  the  estimation  of  the  actual  state  of 
charge  (SOC)  and  state  of  health  (SOH)  of  the  battery  is  rather 
crucial  to  predict  e.g.  the  availability  of  power  and  energy  in 
hybrid  electric  vehicles  (HEV)  and  electric  vehicles  (EV).  During 
lifetime,  the  resistance  as  well  as  the  capacity  and  in  consequence 
SOC  of  every  lithium-ion  cell  changes  through  electrochemical 
degradation  processes  like  electrolyte  decomposition  or  growth  of 
a  solid  electrolyte  interface  (SEI)  on  the  anode  surface  [1],  This 
variation  is  hard  to  measure  directly  in  a  vehicle,  but  has  to  be 
known  for  an  accurate  range  and  power  prediction.  Without  a  SOH 
correction  or  update  by  the  battery  management  system  (BMS), 
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the  driver  will  experiences  an  overestimated  range  or  less 
acceleration. 

Because  of  this  necessity  of  a  reliable  SOC  and  SOH  prediction, 
there  were  a  lot  of  publications  in  the  last  few  years.  The  appli¬ 
cation  of  Kalman  filtering  can  be  regarded  as  a  state  of  the  art 
technique.  Usually  Extended  Kalman  filtering  is  used  (e.g.  Ref.  [2]), 
but  also  Unscented  or  Sigma-Point  Kalman  filtering  is  utilized  (e.g. 
Ref.  [3]). 

The  use  of  machine  learning  techniques,  and  especially  support 
vector  machines,  has  grown  in  popularity  over  the  last  few  years, 
not  only  in  SOC  and  SOH  estimation  [4—6]  but  in  various  applica¬ 
tion  areas.  A  good  overview  of  the  state-of-art  battery  prognostic 
algorithms  is  given  by  Zhang  and  Lee  [7],  where  not  only  Kalman 
filters  and  support  vector  machines,  but  also  other  techniques,  such 
as  neural  networks  or  particle  filters  are  presented. 

This  paper  is  based  on  the  master  thesis  of  Christian  Appel  and 
presents  a  novel  method  for  SOC  and  SOH  prediction  using 
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a  standard  Kalman  filter  for  linear  systems,  an  Unscented  Kalman 
filter  and  machine  learning  techniques  is  presented.  Due  to  its 
general  structure  and  design;  this  approach,  demonstrated  at 
a  nickel  manganese  cobalt  (NMC)  pouch  cell,  is  valid  independent  of 
cell  chemistry  and  design.  Even  if  no  direct  understanding  of  the 
ageing  mechanism  or  sources  can  be  gained  by  such  SOH  algorithm, 
information  can  be  collected  about  the  caused  degradation  of 
particular  ageing  factors.  Thus,  the  main  factors  can  be  identified 
and  minimized  by  an  optimized  operation  strategy  on  the  BMS 
later  on. 

After  the  introduction  and  the  presentation  of  the  state  of  the  art 
SOC  and  SOH  estimation,  a  brief  description  of  the  mathematical 
methods  used  in  this  approach  is  given  in  the  second  section, 
followed  by  the  description  of  the  measurements  and  ageing  test. 
The  fourth  section  describes  the  implementation  of  the  battery 
model  and  the  application  of  the  mathematical  methods  and 
subsequently  the  results  are  presented  and  discussed. 


2.  Mathematical  methods 

In  this  section  a  description  of  the  mathematical  methods 
applied  in  this  study  is  given.  Namely  this  is  a  standard  Kalman 
filter,  an  Unscented  Kalman  filter  and  the  support  vector 
regression. 


2.1.  Kalman  filter  for  linear  systems 

The  first  used  technique  will  be  the  well-known  Kalman 
filter  (KF)  for  linear  systems,  developed  by  Kalman  in  1960  [8], 
A  good  introduction  to  Kalman  filters  can  be  found  in  Ref.  [9], 
Thus,  only  a  presentation  of  the  applied  algorithm  is  given  in 
this  paper. 

A  n-dimensional  linear  dynamic  discretized  system  of  a  discrete 
random  process  consists  of  the  following: 

(i)  A  process  model: 

*k+i  =  Ak*k  +  Bkuk  +  vk,  (1) 

where  xk  is  a  discrete  random  process,  which  is  called  state  vector, 
Ak  e  IRnxn  the  state  transition  matrix,  Bk  e  IRnxm  the  control  input 
matrix,  uk  e  IRm  a  deterministic  control  vector  and  vk  an  uncorre¬ 
lated  discrete  random  process,  which  is  called  process  noise, 
where: 

Qk  =  E(vkvTk) 

(ii)  An  observation  model: 

zk  =  Hkxk  +  wk,  (2) 

where  Hk  e  [R,xn  is  the  observation  matrix  and  wk  an  uncorrelated 
discrete  random  process,  which  is  called  measurement  noise, 
where: 

Rk  =  E(wkwl) 

Presuming  that  the  following  assumptions  hold: 

E(wkvj)  =  0  Vj,  k 


E(wkxJ )  =  0  for  j  <  k 

(4) 

E(vkxj |  m  0  for  j  <  k 

(5) 

Let  xk  =  x(k\.k)  the  a  posteriori  state  estimate  and 
Pk  —  P(k\k)  be  the  a  posteriori  error  covariance  matrix  at  time 
step  k,  given  measurements  up  to  time  step  k,  then  the  standard 
Kalman  filter  equations  for  the  described  linear  system  are  the 
following: 


(i)  Calculate  state  estimate  and  estimate  covariance 
prediction  (a  priori): 

*(k+l|k)  =  Akx(k\k)  +  Bkuk 
P(k+  l|fc)  =  AkP(k\k)A[  +  Qj, 

(ii)  Calculate  innovation  and  innovation  covariance: 

rk+ 1  =  zk+i  -Hk+1x(k+l\k) 

Sk+ 1  =  HknP(k+)\k)HTkn+Rkli 

(iii)  Calculate  Kalman  gain: 

Kk+1  =P(k+l|k)H[+1%Jt 

(iv)  Calculate  state  estimate  and  estimate  covariance 
update  (a  posteriori): 

x(k  +  l|k  +  1)  =  x(k+  l|k)  +  /4+irfc+i 
P(k  +  1  |k  +  1)  =  P(k  +  1 1  k)  +  Kfc+1Sk+1K[+1 


2.2.  Unscented  Kalman  filter 

For  our  purpose  a  nonlinear  extension  to  the  Kalman  filter  is 
applied:  the  Unscented  Kalman  filter  (UKF),  using  the  unscented 
transformation  as  described  in  Ref.  [10].  Rather  than  linearizing  the 
nonlinear  functions  as  in  the  Extended  Kalman  filter  (EKF),  the 
probability  distribution  will  be  approximated  by  a  certain  number 
of  sigma  points  and  those  will  be  transformed  by  the  nonlinear 
system  functions  to  approximate  the  mean  and  covariance 
estimate. 

The  discretized  general  nonlinear  dynamic  system  has  the 
following  form,  analogous  to  the  linear  systems  (1)  and  (2)  and 
requiring  assumptions  (3)— (5)  still  hold: 

(i)  Process  model: 


Xfc+1  =  fk(xk,uk,vk)  (6) 


(ii)  Observation  model: 

(3)  zk  =  hk(xk,wk) 


(7) 
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(i)  Choose  n  +  1  sigma  points  as  follows: 

x(°)  =  X  =  E(x) 

x®  =  x(°)  +  (y/(n  +  k)P)J,  i  =  1 

x(n+i)  =  x(0)  _  (^(n  +  K)p)  *  i  =  1, 

where  (\/(n  +  k)P),  is  the  ith  row  of  the  matrix  square  root 
(x/(n  +  x)P). 

Weights  of  the  sigma  points  are  chosen  as  follows: 


e  estimate  prediction  (a  priori): 


X®(k+Hk)  =  f(x®(k\k),Uk)+vk 
2n 

x(k  +  l|k)  =  J2  W(i)x(i)(k+ l|fc) 

i  =  0 


(iii)  Calculate  estimate  covariance  prediction  (a  priori): 

P(k  +  1 1 k)  =  0*  +  W®  |x®(k  +  l|fc)  -  x(k  +  1  |k)J 

x  [x®(fc+l|fc)-x(k  +  l|k)j 

(iv)  Calculate  innovation: 


Z®(fc  +  l|fc)  =  h(x®(k  +  \  \k),uk)  +wk 

z(k  +  l|k)  =  £  W®z©(k  +  l|k) 

1  =  0 


(v)  Calculate  Kalman  gain: 


Pyy(k  +  1  \k)  =  Rk  +  f)  If/®  \z®(k  +  1 1 k)  -  z(k  +  1  |k)l 
(k  +  1  |k)  —  z(k  +  ll^  Pxy(k+  l|k) 

=  W®  |x®  (k  + 1  |k)  —  x(k  + 1  |k) j  |z®  (fc  +  l[fc) 

—  z(k+  l|k)J  T/Cfc 

=  Pxy(fc  +  1  |/<)Pyy(fc  +  life)1 

(vi)  Calculate  state  estimate  update  (a  posteriori): 

x(k  +  l|k  + 1)  =  x(k  +  l|k)  +  Kk(zk  -  z(k+  l|k)^ 

(vii)  Calculate  estimate  covariance  update  (a  posteriori): 
P(k  + 1  |k  + 1)  =  P(k+l|k)  +  /fkPyj,(k+l|k)/Cfcr 


23.  Support  vector  regression 

The  second  method  used  in  this  approach  is  the  support  vector 
regression  (SVR).  A  comprehensive  introduction  can  be  found  in 
Ref.  [11],  Given  a  set  of  training  data: 

{(xi,yi),(x2,y2),...,(x„,yn)}cRn  x  r  (8) 

where  x,-  e  Rn  are  the  input  parameter  vectors  and  y,-  e  K  are  the 
target  values,  the  task  in  the  so-called  r-incentive  support  vector 
regression  is  to  find  a  linear  function 

f(x)  =  (w,x)  +  b  (9) 

as  flat  as  possible  in  the  way  that  the  values  ^x,)  have  at  most  e 
deviation  from  the  targets  y,.  Introducing  “slack  variables"  f;,  ?,■  to 
create  a  “soft  margin”  and  therefore  allowing  measurement  errors 
and  be  able  to  cope  with  otherwise  infeasible  constraints,  leads  to 
a  dual  optimization  problem.  Utilizing  the  Karush— Kuhn— Tucker 
(KKT)  condition,  the  parameter  b  in  function  f[x)  can  be  calcu¬ 
lated.  Obviously  SOC  and  SOH  prediction  will  not  be  a  linear 
regression  problem,  so  a  method  for  the  nonlinear  case  is  required. 
Hereby,  the  idea  is  to  transform  the  input  data  nonlinearly  in 
a  higher  dimensional  feature  space  F : 

<b  ■  IRn— »F 

x^4>(x) 

The  equations  in  the  dual  optimization  problem  are  solely 
dependent  on  inner  products  of  the  form  (X;,Xj).  Hence  in  feature 
space  it  is  sufficient  to  know 

k(Xj,Xj)  =  {<j>{Xi),(j>(Xj))  (10) 

with  a  so  called  kernel  function  k  to  represent  an  inner  product  in 
feature  space  and  the  nonlinear  transformation  0  has  not  be  known 
explicitly.  This  is  the  kernel  trick  and  Mercer’s  theorem  (1909)  gives 
a  condition  for  such  functions. 

In  this  study  an  E-incentive  support  vector  regression  is  applied 
and  Gaussian  Kernels  of  the  following  form  will  be  used: 

(11) 

3.  Experiments 

In  this  section  a  description  of  the  measurements  taken 
for  the  present  study  and  details  of  the  used  cell  type  are 
given. 

High-energy  lithium-ion  pouch  cells  with  a  graphite  anode  and 
a  NMC  cathode  were  used.  The  cells  have  a  nominal  voltage  of  3.6  V 
and  a  nominal  capacity  of  10  Ah. 

In  order  to  have  the  chance  to  validate  the  subsequent  esti¬ 
mation,  all  cells  were  characterized  by  electrical  parameter  tests. 
Besides  a  capacity  test  with  a  C-rate  of  one,  a  pulse  power 
characterization  profile  (PPCP),  based  on  the  “Hybrid  Pulse 
Power  Characterisation  Test”  of  PNGV/FreedomCar  and  the 
EUCAR  “Open  Circuit  Voltage  and  Power  Determination  Test” 
[12],  is  part  of  this  characterization.  It  consists  of  a  certain  pulse 
power  profile  at  different  SOC  steps,  from  90%SOC  down  to  10% 
SOC  in  10%SOC-steps.  The  pulse  power  profile  procedure  can  be 
seen  in  Table  1. 

The  SOC-OCV  curves,  which  are  later  required  in  the  process 
model  of  filter  2,  are  obtained  from  the  PPCP  measurements  as  well 
as  a  time-dependent  resistance 
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which  can  be  calculated  for  each  SOC  step. 

The  recorded  measurement  data  can  be  divided  into  two 
different  classes: 

(i)  Long-term  cell  ageing  tests:A  large  amount  of  storage  and 
cycle  ageing  tests  were  performed  over  more  than  one  year 
with  specific  lithium  ion  cells.  During  the  test  duration,  regular 
check-ups  were  conducted  to  record  the  capacity  fade  in  the 
cells.  The  check-ups  consist  of  standardized  capacity  and 
power  tests  as  described  above. 

(ii)  Short-term  driving  cycle  cell  measurements :The  cells,  which 
passed  the  ageing  tests  (i)  and  consequently  have  a  different 
SOH,  were  used  for  two  types  of  driving  cycle  measurements 
based  on  the  FTP-72  cycle  to  simulate  real  world 
requirements: 

(ii.a)  Full-charged  cells  are  completely  discharged  by  FTP-72 
load  profiles. 

(ii.b)  Full-charged  cells  are  discharged  down  to  ca.  40%  by  a  load 
profile  consisting  of  FTP-72  cycles,  pauses  and  a  charge 
period. 

The  load  profile  based  on  the  FTP-72  cycle  and  the  load  profile  of 
the  combined  cycle  measurements  can  be  seen  in  Figs.  1  and  2, 
respectively.  While,  the  first  type  of  measurements  is  only  used  for 
having  cells  at  different  SOH  with  a  known  history  and  for  the  SVR 
application.  The  second  type  is  used  for  a  real-time  SOC  and  SOH 
estimation. 


Fig.  2.  Applied  current  load  profile  for  combined  cydic  measurements  consisting  of 
pauses,  drive  cydes  and  recharging  with  a  total  time  of  7  h. 


where  Uocv.fc  is  the  open-current  voltage  (OCV)  at  time  step  4-  The 
OCV  corresponds  to  a  certain  SOC,  which  is  implemented  by  a  non¬ 
linear  characteristic  curve  Uocv(SOQj).  Moreover,  this  curve  is 
assumed  to  be  consistent  over  lifetime,  if  the  SOC  calculation  is 
based  on  the  actual,  not  the  nominal,  capacity.  Furthermore  the 
over-voltages  Ui  and  U2  can  be  expressed  by: 


where  Atk+i:  =  4+ 1  -  4  is  the  time  span  between  time  step  k  +  1 
and  k. 

These  equations  will  be  implemented  in  both,  the  standard 
Kalman  filter  (filter  1)  and  the  Unscented  Kalman  filter  (filter  2). 
The  process  model  equation  for  the  state  of  charge  is: 


4.  Implementation 

In  this  section  the  implemented  battery  model  is  introduced  and 
the  derivation  of  the  Dual  Kalman  filter  (DKF)  and  the  application  of 
the  SVR  algorithm  is  presented. 

The  underlying  battery  model  is  a  state  of  the  art  electric 
equivalent  circuit  model  with  2  RC  members  and  is  shown  in 
Fig.  3  [13], 


Fig.  1.  Applied  current  load  profile  based  on  FTP-72  drive  cycle. 


SOCfc+1  =  SOCfc  +  /k'Atkl1fOO  (15) 

Lcell 

where  4  is  the  measured  current  at  time  step  k,  which  is  constant 
on  the  time  interval  [4,4+1]. 

In  this  study,  a  combination  of  a  standard  Kalman  filter  and  an 
Unscented  Kalman  filter  is  used.  The  first  filter  estimates  the 
polarization  over-voltage  U 1,  the  diffusion  over-voltage  U2  and  the 
ohmic  resistance  R0hm-  The  output  of  this  filter  is  fed  into  the  second 
filter,  estimating  SOC  as  well  as  the  polarization  and  diffusion 
resistances  Ri  and  R2.  In  the  system  model  of  the  second  filter,  the 
actual  cell  capacity  estimate  of  the  SVR  can  be  used  to  predict  the 
SOC. 


Ri  R  2 


Ui  Ui 


Fig.  3.  Used  electric  equivalent  circuit  model  of  a  lithium-ion  cell  and  corresponding 
internal  cell  parameters. 


Table  2 

Equations  for  filter  1  (KF). 
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Control  input 
Observation 


Uk  =  4 

Zk  =  UceU,k  -  Uocv,k 


Process  model  xk+1  =  Ak-xk  +  Bk  uk 

Observation  model  zk  =  (1  1  Ik)-Xk 


The  motivation  for  this  Dual  Kalman  filter  is  based  on  three 
advantages: 

•  Decoupling  of  estimations,  thus  reduction  of  interactions  and 
avoiding  of  building-ups  of  a  filter 

•  Separation  of  variables,  which  cannot  be  estimated  by  a  single 
filter 

•  Reduction  of  computation  efforts,  since  two  filters  of  lower 
dimensions  are  faster  than  one  higher  dimensional  one  [14] 


£fc  - >  tk+ 1 


Both  filters  use  the  input  signal  current  and  the  measured  cell 
voltage.  Measurement  data  were  obtained  by  the  previously 
described  driving  cycle  tests  to  simulate  real  world  requirements. 
Therefore  the  dual  filter  can  be  used  online,  that  means  e.g.  it  can  be 
implemented  in  the  battery  management  system  to  predict  the 
power  and  energy  availability. 

The  equations  for  both  filters  can  be  found  in  Tables  2  and  3, 
respectively.  The  time  constants  x\  =  RiCi  and  12  =  R2C2  were 
assumed  to  be  constant  in  both  filters.  This  assumptions  was  rather 
made  than  keeping  the  capacities  Ci  and  C2  constant,  because  of  the 
fitting  results  and  robustness  considerations.  The  dual  filter  layout 
can  be  found  in  Fig.  4. 

The  actual  cell  capacity  is  also  required  in  the  process  model  of 
the  second  (Unscented)  Kalman  filter.  It  is  possible  to  take  the 
nominal  capacity  or  any  other  recent  measurement  result.  In  this 
study,  the  actual  relative  cell  capacity  Crei,fc  =  Cceiu/CBOT  is  esti¬ 
mated  by  c-SVR,  where  Cbot  is  the  begin-of-test  (BOT)  capacity,  and 
serves  as  an  input  value  for  the  model  parameter  in  the  second 
Kalman  filter.  The  data  for  the  SVR  was  obtained  from  the  long¬ 
term  cell  ageing  tests  (i). 

The  composition  of  the  input  vector  can  be  seen  in  Table  4  with 
an  additional  scaling  to  [-1,1]  for  temperatures  and  [0,1]  for  the 
other  features  [15]. 

Multiple  data  vectors  were  created  for  every  check-up  of  every 
cell,  such  that  every  time  interval  between  a  check-up  and  all 
previous  check-ups  including  the  BOT  measurement  is  accounted. 


Equations  for  filter  2  (UKF). 


In  this  way  a  number  of  1899  training  vectors  and  566  testing 
vectors  was  obtained. 

The  choice  of  the  required  e-SVR  parameters  can  be  seen  in 
Table  4,  where  the  value  of  e  was  chosen  because  of  the  expected 
measurement  error  and  cell  spread.  The  kernel  function  parameter 
cr  =  1  /2<r2  and  the  regularization  constant  C  were  obtained  by  grid 
search  and  5-fold  cross  validation  [15]  based  on  the  training  data. 

The  SVR  part  is  used  offline  and  the  ideas  are  that  new  training 
data  can  be  fed  into  the  model  e.g.  regularly  when  a  full  charge  is 
performed  or  at  regular  services. 

In  addition,  the  SVR  could  also  be  used  to  obtain  initial  values  for 
the  resistances  estimated  in  the  DI<F  or  to  deliver  estimates  of  the 
state  of  health  on  its  own. 

A  very  interesting  further  application  of  the  SVR  would  be  to 
determine  the  parameters  of  an  equation  describing  the  ageing  or 
the  lifetime  instead  of  predicting  the  capacity  directly.  A  possible 
equation  would  be  e.g: 


^actual  —  Cbol-  (lOO  a-fl3  t? '^throughput)  (16) 


Table  4 


Data  structure  and  set  parameters  of  the 

;  support  vector  regression. 

Variable 

Description 

(a)  Composition  of  the  data  vectors 
Input  vector 

Crei.fi 

Tmin 

Relative  cell  capacity  at  time  ti 

Tmax 

Maximal  temperature 

d  =  t2-t, 

Time  span 

Q 

Mean  current  rate 

Amount  of  charge  and  discharge 

SOCmax 

Maximal  SOC 

Asoc 

Target  value 

SOClift 

CrelA 

Relative  cell  capacity  at  time  t2 

SVR  parameter 

(b)  Initialised  parameters 

Kernel  function  k(x,y) 

where  112a2  =  2.9 

£ 

0.005 

c 

60 

Zfc  = 
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to  describe  the  capacity  fade  over  calendar  and  cycle  life.  The 
dependency  of  the  parameters  on  the  operation  conditions  could 
now  be  obtained  by  the  SVR.  This  approach  was  tested  and  showed 
also  very  satisfying  results  for  the  training  as  well  as  the  validation 
with  unknown  ageing  conditions.  Two  main  advantages  can  be 
found  for  this  indirect  approach.  First,  measurements  inaccuracies 
are  compensated  and  the  estimation  is  more  robust  against 
measurement  errors  or  noise.  Second,  by  obtaining  these  parame¬ 
ters  an  additional  prognosis  is  possible.  Especially,  if  a  large  number 


of  ageing  test  is  available,  the  parameters  are  well  suited  for  a  life¬ 
time  prediction. 

In  this  study,  the  SOH  is  not  predicted  directly,  but  can  be 
calculated  from  the  estimated  actual  cell  capacity  and  the  ohmic 
resistance  defining  specific  end-of-life  criteria.  Here,  the  SOH  will 
be  defined  the  following  way: 

!1 00%  begin  -  of  -  life(BOL) 

0%  end  -  of  -  life(EOL) 
linear  interpolated  in  between 

with  the  following  typical  automotive  EOL  criteria: 


Fig.  5.  During  a  drive  cycle  estimated  ohmic  resistance  and  over  voltages  of  four  cells 
with  different  capacity  SOH  (blue  line  94.5%,  green  line  90.0%,  red  line  95.2%,  light  blue 
line  98.0%)  by  the  first  Kalman  filter.  (For  interpretation  of  the  references  to  colour  in 
this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Reol  =  2-Rbol 
Qol  =  0.8-Cbol 


Fig.  6.  During  a  drive  cycle  estimated  SOC,  polarization  resistance  and  diffusion 
resistance  of  four  cells  with  different  capacity  SOH  by  the  second  Kalman  filter. 
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Both,  the  Dual  Filter  and  the  support  vector  regression  have 
been  implemented  in  MATLAB  code,  where  a  common  SVM 
toolbox,  LIBSVM  [16],  is  used  for  the  SVR  part. 


5.  Results  and  discussion 

The  driving  cycle  measurement  data  (ii.a)  is  first  applied  to  the 
dual  filter  algorithm.  The  estimated  states  over  time  for  four 
different  cells  are  plotted  in  Figs.  5  and  6  for  filter  1  (KF)  and  filter  2 
(UKF),  respectively.  It  can  be  seen,  that  the  ohmic  resistance  esti¬ 
mate  Rohm  varies  clearly  between  the  differently  aged  cells  and  it 
increases  as  expected  with  decreasing  state-of-charge. 

The  SOC  estimate  shows  very  uniform  characteristics  starting 
with  a  value  100%SOC  and  ending  with  almost  0%SOC.  Especially, 
comparing  the  curve  of  cell  two  and  three,  a  different  cell  capacity 
can  be  guessed. 

Also,  the  evolution  of  the  over-voltages  Ut  and  U2  over  time 
seem  feasible,  but  cannot  be  validated  directly.  Flowever,  two 
interesting  points  are  obviously.  First,  the  same  trends  of  the  over¬ 
voltages  is  found  for  all  four  cells  and  second,  the  over-voltage  due 
to  diffusion  is  four  times  higher  than  due  the  polarization  effect.  For 
the  resistances  R\  and  R2  a  little  noisy  behaviour  in  a  fixed  constant 
range  of  0.2-0.3  mQ  and  4-5  mQ,  respectively,  independently  of 
initial  values  or  the  type  of  driving  cycle  measurement  is  found.  In 
contrast  to  Rohm,  no  significant  influence  of  the  SOFI  is  visible. 

Since  the  concrete  values  of  the  over-voltages  U\  and  U2  and  the 
corresponding  resistances  Ri  and  R2  are  not  essential,  only  the  SOC 
and  Rohm  estimates  for  the  second  driving  cycle  measurement  (ii.b) 
are  shown.  The  evolution  over  time  for  those  two  states  is  given  by 
Fig.  7a  and  b,  respectively. 

Again,  the  ohmic  resistance  of  the  investigated  cells  shows 
similar  but  shifted  curves.  Also  the  SOC  of  the  four  cells  over  the 
drive  cycle,  starting  at  different  times  due  to  visibility,  show  almost 
identical  trends.  However,  a  deeper  look  at  the  reached  SOC  level 
during  both  pauses  shows  slightly  different  values.  This  is  caused 
by  the  different  capacity  SOH  of  the  cells  and  can  be  evaluated.  To 


Fig.  8.  Comparison  of  calculated  (dashed  line  with  circles)  and  estimated  (solid  lines 
with  triangles)  ohmic  resistance  of  three  cells  with  different  resistance  SOH  (red  line 
127%,  green  line  108%,  blue  line  101%).  (For  interpretation  of  the  references  to  colour  in 
this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


validate  the  results  of  applied  dual  filter,  the  PPCP  test  of  Section  3  is 
used  including  the  SOC— OCV  curves  as  well  as  the  time  dependent 
resistance  of  Eq.  (12)  of  each  cell.  Since  the  voltage  drop  caused  by 
the  ohmic  resistance  happens  instantaneously,  one  tries  to  account 
the  voltage  difference  of  the  smallest  analysable  time  difference.  In 
this  case,  the  smallest  time  residual  for  the  measurement  was  0.5  s 
and  will  contain  also  polarization  effects.  Nevertheless,  this  resis¬ 
tance  is  in  a  small  range  as  demonstrated  in  Fig.  6  and  consequently 
R0.5  s  provides  a  good  approximation  for  the  ohmic  resistance. 

It  is  therefore  possible  to  validate  the  state  estimates  SOC  and 
Rohm  with  the  measured  open-current  voltages  after  relaxation  and 
the  calculated  R0.5  s  resistances,  respectively.  Obviously,  only  the 
SOC  estimates  at  the  rest  periods  of  the  second  drive  cycle 
measurement  (ii.b)  can  be  validated.  The  ohmic  resistance  estimate 
Rohm  will  be  validated  comparing  the  results  of  the  first  driving 
cycle  measurement  (ii.a),  because  it  covers  the  whole  SOC  range 
from  100%SOC  to  almost  0%SOC  with  the  calculated  R0  5  s  resistance 
based  on  the  PPCP  measurements. 

A  comparison  of  the  calculated  and  the  estimated  resistance 
values  over  state  of  charge  for  three  selected  cells  is  given  by  Fig.  8. 
Both  resistances  show  similar  characteristics  with  a  slightly 
different  gradient.  The  calculated  resistances  show  a  higher 
increase  with  decreasing  SOC  compared  to  the  estimates.  This  is 
explainable  with  the  voltage  drop  caused  by  the  OCV  decrease 
when  discharging  the  cell,  which  is  accounted  in  the  dual  filter 
model  but  not  in  the  calculation.  The  lower  the  SOC,  the  higher  the 
OCV  decrease  and  therefore  the  difference  in  the  resistances. 


Table  5 

Estimates  SOQ,a  and  deviation  ASOC  from  the  measurement  value  SOCactuai  for 
different  start  values  of  Ccell. 


Cell  1  Cell  2  Cell  3 


Cell  4 


(a)  Estimated  SOC  and  deviation  ASOC  from  the  measured  value . 


of  the  first  rest  period 
SOCact 

SOCest(Qell,a) 


SOCest(Ccell,l  =  8  Ah) 
ASOC 

SOCeflfCcdu  =  12  Ah) 
ASOC 

(b)  Estimated  SOC  and  1 


7235  70.85 

71.79  70.33 

0.56  0.52 

71.64  70.20 

0.71  0.65 

71.91  70.49 

0.44  0.36 

ition  ASOC  from  the 


72.96 

72.66 

0.30 

72.47 

0.49 

72.78 

0.18 


measured  value . 


of  the  second 
SOCact 

SOCest(Ccelu) 

ASOC 

SOCestCQeUj  =  8  Ah) 
ASOC 

SOCes,(Cceiii2  =  12  Ah) 
ASOC 


65.06  63.06 

65.04  62.88 

0.02  0.18 

64.84  62.70 

0.22  0.36 

65.19  63.10 

0.13  0.04 


65.80 

65.70 

0.10 

65.49 

0.31 

65.84 

0.04 


73.55 

72.86 

0.69 

72.60 

0.95 

72.97 

0.58 

:  the  end 


65.40 

65.35 

0.05 

65.13 

0.27 

65.45 

0.05 
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Data  sets 


Data  sets 

Fig.  9.  Validation  of  the  SVR  training  (a)  and  capacity  prediction  (b)  by  comparison  of 
SVR  (blue  line)  to  calendar  and  cycling  ageing  measurements  (green  line),  (a)  Training 
of  the  SVR.  (b)  Prediction  by  SVR.  (For  interpretation  of  the  references  to  colour  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Subsequently,  a  very  good  quality  of  the  resistance  estimate  can  be 
followed  too. 

Table  5  shows  the  SOC  estimates  for  three  different  values  for 
the  cell  capacity  model  parameter  CCeii,  where  CCeii,a  is  the  actual 
“true”  cell  capacity  obtained  from  capacity  tests  and  Cceii,i  and  CCeii,2 
are  fixed  values,  which  are  up  to  20%  lower  and  higher,  respectively. 

It  can  be  seen,  that  the  maximal  deviation  of  the  SOC  estimate 
from  the  true  value  is  under  1%,  even  if  the  model  parameter  cell 
capacity  Cceii  has  20%  deviation  from  the  true  value.  By  assuming 
again  a  consistent  form  of  the  OCV-curve  over  lifetime,  the  actual 
cell  capacity  can  now  be  calculated  easily  by  using  the  OCV/SOC 
relation  and  the  discharged  Ah-amount. 

Now,  the  possibility  to  obtain  an  initial  Cceii  for  the  DKF  by  the 
use  of  SVR,  is  tested  and  validated  with  regard  of  a  maximum  error 
limit  of  20%  for  an  unknown  cell  history.  As  mentioned  before, 
a  total  of  almost  two  thousand  data  vectors  from  ageing  tests  (i) 
was  applied  in  order  to  train  the  SVR.  A  very  successful  training  is 
achieved  as  evidenced  by  Fig.  9a.  The  curve  of  all  storage  and 
cycling  ageing  tests  is  reproduced  by  the  SVR  throughout  with 
a  high  accuracy.  Afterwards,  the  SVR  model  is  used  to  predict  six 


ageing  scenarios  which  were  not  part  of  the  training.  Also  this 
validation  shows  a  high  accuracy  as  illustrated  in  Fig.  9b.  The  overall 
trend  of  all  cells  is  reproduced  correctly  and  just  small  deviations, 
with  an  average  error  in  the  range  of  the  cell  spread  between  the 
three  tested  cells,  can  be  found.  Looking  at  the  maximum  error, 
even  for  the  strongest  degradation  a  value  below  20%  is  reached. 
Thus,  the  predicted  capacity  is  outstandingly  suitable  to  be  used  as 
an  input  value  for  the  dual  Kalman  Filter. 

6.  Conclusions 


In  this  study  a  robust  and  powerful  real-time  SOC  and  SOH 
estimation  for  lithium-ion  batteries  has  been  developed.  It  was 
shown,  how  the  mathematical  methods,  namely  minimum  vari¬ 
ance  estimation  and  machine  learning,  can  be  combined  for  this 
purpose.  The  performance  of  the  resistance  and  SOC  estimation 
have  been  validated  using  results  of  pulse  power  characterization 
tests  and  the  open-current  voltage  curves,  respectively. 

The  accuracy  of  the  SOC  and  resistance  estimation  can  be 
regarded  superior  to  existing  publications.  Even  with  a  capacity 
error  of  20%,  the  DKF  showed  excellent  results  with  an  SOC  esti¬ 
mation  error  below  1%  as  well  as  comprehensible  values  for  all 
resistances  in  case  of  realistic  drive  cycles.  Moreover,  it  was  showed 
that  the  SVR  can  be  used  offline  to  gain  a  convenient  initial  capacity 
for  the  DKF  or  even  for  lifetime  prediction. 

In  a  next  step,  the  realization  of  the  algorithms  on  a  BMS  should 
be  evaluated  in  terms  of  computation  and  accuracy  for  vehicle  data 
instead  of  laboratory  data. 
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